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Abstract 

A typical gap noise problem is considered in 
this paper using the new space-time conser- 
vation element and solution element (CE/SE) 
method. Implementation of the computation is 
straightforward. No turbulence model, LES 
(large eddy simulation) or a preset boundary 
layer profile is used, yet the computed fre- 
quency agrees well with the experimental one. 

1 Introduction 

Gap noise is an important issue in aerospace and auto 
industries. For example, modem automobiles are re- 
quired to be comfortable and quiet, however, due to vor- 
tex shedding in high speed air flow over door gaps, un- 
pleasant tonal noise may be generated. Strong oscilla- 
tions occur in a feedback cycle in which the vortices 
shed from the upstream edge of the gap convect down- 
stream and impinge on the other edge, generating acous- 
tic waves that propagate upstream to excite new vortices. 
Numerical simulation of such a complicated process re- 
quires a scheme that can : (a) resolve acoustic waves 
with low dispersion and numerical dissipation, (b) han- 
dle nonlinear and discontinuous waves (e.g. shocks), and 
(c) have an effective (near field) non-reflecting boundary 
condition (NRBC). The new space time conservation el- 
ement and solution element method, or CE/SE for short, 
is a numerical method that meets the above requirements 
[1-4]. A detailed description of the 2-D CE/SE Euler 
scheme can be found in [ 1 , 2 ], only a brief sketch is given 
here. 

In nature, the new method may be categorized as a 
finite volume method, where the conservation element 
( CE ) is equivalent to a finite control volume (or cell) 
and the solution element ( SE ) can be understood as the 
cell interface. However, due to its rigorous treatment of 
the fluxes and geometry, it is different from the existing 
schemes. As demonstrated in the previous papers [4-7], 
the CE/SE scheme features: 
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( 1 ) space and time treated on the same footing, the in- 
tegral equations of conservation laws are solved for with 
second order accuracy, 

( 2 ) high resolution, low dispersion and low dissipa- 
tion, 

( 3 ) novel, truly multi-dimensional, simple but effec- 
tive non-reflecting boundary condition, 

(4) Straight forward implementation of computation, 
no numerical fix or parameter choice is needed; 

( 5 ) robust enough to cover a wide spectrum of com- 
pressible flow: from weak linear acoustic waves to 

strong, discontinuous waves (shocks), appropriate for 
linear and non-linear aeroacoustics. 

Currently, the CE/SE scheme has been developed to 
such a stage that a 3-D unstructured CE/SE Navier- 
Stokes solver is already available. However, in the 
present paper, as a general introduction to the CE/SE 
method, only the 2-D unstructured Euler CE/SE solver 
is chosen as a prototype and is sketched in Section 2. 
Then application of the CE/SE scheme to the auto-door 
gap noise problem is depicted in Sections 3. 

2 The 2-D Unstructured CE/SE Scheme 

The CE/SE scheme can be used with either structured 
or unstructured grids. Here, the unstructured CE/SE is 
used as a prototype for introductory description. 

2.1 Conservation Form of the Unsteady Euler 
Equations 

Consider a dimensionless conservation form of the un- 
steady Euler equations of a perfect gas. Let p, u, v y 
p, and 7 be the density, streamwise velocity compo- 
nent, transversal velocity component, static pressure, and 
constant specific heat ratio, respectively. The 2-D Eu- 
ler equations then can be written in the following vector 
form: 

Ut+Fs-f-G^O, (1) 

where x, y, and t are the streamwise and transversal co- 
ordinates and time, respectively. The conservative flow 
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variable vector U and the flux vectors F and G, are given 

by: 


U = 




G = 


( Gl \ 

G 2 

G a 


\gJ 


a node ( e.g., D,E,F in Fig. 1 ) is here a set of interface 
planes in £3 that passes through this node (e.g. planes 
DAA f D ( , DBB'D DBOA associated with node D). 
Within a given solution element SE(j, n), where j, n 
are the node index,and time step respectively, the flow 
variables are not only considered continuous but are also 
approximated by linear Taylor expansions: 


with 

r r i = p , U 2 - pu , U 3 = pi\ U\ = pf ( 7 — 1 )-t-p(u 2 -M; 2 )/ 2 ; 

*i =r 2 , F^Lhlh/Uu 
r -2 = (7 - 1)^.1 + [(3 - 7 )ui - (7 - i)f/|] /Wl 
Fa = 1U2U4/U1 - (0 - 1 ) 0 2 [ 0 * + r/f] / 2 C/, 2 , 

Gl = [ 3* G2 — r r 2 £ 3 / £ r i , 

Ga = (7 - l)t'4 + [(3 - 7)C'I - (7 - 1 )Ul] /2U U 

G 4 = yUaUt/Ui - (7 - 1)03 [0? + Ul] /2Uf. 

By considering (x>y,t) as coordinates of a three- 
dimensional Euclidean space £3 and using Gauss’ diver- 
gence theorem, it follows that Eq. (1) is equivalent to the 
following integral equation: 

<£ H m dS = 0 , m = 1,2,3, 4, (2) 

JS(V) 

where S(V) denotes the surface around a volume V in 

E-s and H m = (F ra ,G m ,l7 m ),m = 1,2, 3, 4. 

2.2 Unstructured CE/SE 

The CE/SE scheme is naturally adapted to unstructured 
triangle grids. Each triangle center (O in AABC ) 
and its 3 neighboring triangle centers (J9, E, F ) form 3 
conservation elements or CEs (quadrilateral cylinders - 
ADBO, BECO and CFAO), as shown in Fig. 1 . These 
triangle centers are the nodes where the unknowns U, 
U x , U y are defined. 

Assume that at the previous time level n (Fig. 1), U, 
and its spatial derivatives Uj., U y are given at all the 
nodes (triangle centers), the CE/SE time marching is 
based on a ’tripod’ mode, i.e U, U x , at the new 
time level n + 1 (e.g., the shaded circle O f in Fig. 1) are 
computed from their data at its surrounding neighboring 
nodes at time level n (e.g., D in Fig. 1, note that O itself 
is not used). 

During the time marching, the above flux conserva- 
tion relation ( 2 ) in space-time is the only mechanism that 
transfers information between node points. A conserva- 
tion element CE , or computational cell, is the finite vol- 
ume to which (2) is applied. Discontinuities are allowed 
to occur in the interior of a conservation element and 
nowhere else. A solution element SE associated with 


U*(x, y,t;j,n) = TJ] + ( U*)’ l (x - xj)+ 

(U (3) 

F *(*,», = Fj + (F x )”(a- - Xj)~ 

(F„ )''(.</ - Vj) + (F t yj(t - t n ), (4) 

G*( x,y,t;j,n) = G" + (G X )](x - Xj)+ 

(G y )?(y - yj) + (G t yj(t - t% ( 5 ) 

where j is the node index of D.E or F, the partial 
derivatives of F and G can be related to the correspond- 
ing one of U by using the chain rule and U t can be 
obtained from ( 1 ). These Taylor expansions are used 
to accurately evaluating fluxes on SE f s e.g., on planes 
DAA'D\ DBB'D ' or DBOA. 

In principle, the number of equations derived from 
these flux conservation laws (each CE provides 4 scalar 
equations and the total number of scalar equations is 
12 ) matches the number of unknowns (here 12 scalar 
unknowns). All the unknowns are solved for based on 
these relations. No extrapolations (interpolations) across 
a stencil of cells are needed or allowed. 

In practice, there is no need to solve a 12 x 12 equa- 
tion system. In the space-time E% space, (2) is applied to 
the hexagon cylinder ADBECF - A' D 1 B' E'C f F f . U 
at the hexagon center at the new time level is first eval- 
uated and then , U at the center O' of triangle A* B'C ' 
is obtained by Taylor expansion. Since the evaluation 
of U x and involves application of artificial damping 
and some form of limiters (weighted average), the reader 
is referred to the original papers [ 1 -3] of Chang et al for 
details. 

2.3 Non-Reflecting Boundary Conditions 

In the CE/SE scheme, non-reflecting boundary condi- 
tions (NRBC) are constructed so as to allow fluxes from 
the interior domain to smoothly exit to the exterior of 
the domain. Some variants of the NRBC frequently em- 
ployed in CE/SE schemes are: 

(1) Type I - ’steady NRBC’: 

For a ghost grid node ( j, n) lying outside the domain 
at the top (or bottom) of the domain the NRBC requires 
that 

(u x )" = (U„)? = 0, 
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hexagon cytiner ADBECFA - A’D’B’E'C'PA’ and 
Its 3 CEs, OADB-O'A’D’B’, OBEC-O’B’E’C’, 

OCFA - O’C’F’A’ 

Figure 1: CE/SE unstructured grid, showing 3 CEs and 
the hexagon cylinder in E 3 . 

while Uj is kept fixed at the initially given steady bound- 
ary value. 

(2) Type II - ’outflow NRBC’: 

At the downstream boundary, where there are substan- 
tial gradients in y direction, the NRBC requires that 

(u*)7 = 0, uy = u;r 1/2 (u„)y = (u,,)",- 1/2 , 

where f is the index of an interior node closest to the 
boundary ghost node j and U™ and (U y )j are now de- 
fined by simple extrapolation from the interior. This 
NRBC is valid for either supersonic or subsonic flows. 
It should be noted that although these NRBC’s bear sim- 
ilarity to those used in finite difference schemes, the role 
they play is very different. It can be shown that the above 
NRBC’s allow fluxes to smoothly exit to the exterior of 
the domain. 

3 Numerical results 

Figure 2 illustrates the geometric configuration of 
the door gap noise problem, which represents a typi- 
cal aeroacoustic feedback system, producing sustainable 


ioise 



A=1 5.9mm, B=28mm, C=8.76mm, 
D=3.3mm 


Figure 2: auto-door gap problem and the gap dimensions. 

feedback cycle oscillations as a result of vortex shed- 
ding. The problem is given as a benchmark problem 
at the 3rd CAA (Computational Aeroacoustics) Work- 
shop (NASA Glenn , Cleveland, OH, November, 1999) 
and the corresponding experimental data will be released 
later. The actual dimensions of the auto-door gap cavity 
are also prescribed in Fig. 2. The length scale is chosen 
to be twice the dimension C in Fig. 2, and the speed of 
sound and density of the ambient flow are respectively 
the scales for velocity and density. In the current compu- 
tation, there are 48184 triangle elements for the unstruc- 
tured mesh in the computational domain. These trian- 
gles are actually obtained by dividing a rectangular struc- 
tured mesh cell into 4 pieces. The rectangular cell keeps 
a uniform size of Ax = 0.00625 and A y = 0.01136 
around the area of the gap and the interior of the cavity, 
but grows larger near the outer boundaries. At the gap 
and cavity solid walls, no slip boundary conditions are 
applied. 

3.1 High speed case (v = 50.9m/s) 

The mean flow follows the x-direction with Mach num- 
ber M = 0.1497, corresponding to the car speed of 
50.9m/s. Initially, ambient flow condition is imposed 
on the entire flow field. At the inflow, the M = 0.1497 
flow is imposed. At the top and bottom the Type I NRBC 
is used while at the outflow, the Type II NRBC is speci- 
fied. The domain shown in Fig. 3 is exactly the the com- 
putational domain, no buffer zone is used but still, the 
CE/SE NRBC works well. Fig. 3 is taken at time steps 
4000, 14000, and 24000 (with At = 0.00125), it shows 
isobar snapshots of the computed flow field at different 
stages. The generation of vortices and the generation of 
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Figure 3: Isobar snapshots at different time steps. 


(nonlinear) acoustic waves by vortex- gap edge impinge- 
ment, both inside and outside of the cavity, are clearly 
displayed. 

The computation is totally carried out for 225,000 
time steps. Starting from time step 60,000, when the un- 
steady flow is considered as fully developed, the pressure 
history at the mid-point of the left wall of the cavity is 
recorded and provided later for FFT (fast Fourier trans- 
form) analysis. 

Fig. 4 displays the PSD (power spectrum density) in 
log. scale for the time series. The x — axis denotes the 
reduced frequency - the Strouhal numberS*. St = 1 
is equivalent to a frequency of 19,406 Hz. It is ob- 
served that there are several different tone spikes extend- 
ing to very high frequencies. The frequency at the lowest 
tone is about 1 , 839 Hz. Despite the coarse resolution in 
frequency domain, this agrees well with the experimen- 
tal frequency data 1 , 824 Hz provided by the 3rd C AA 
Workshop [8]. 

Fig. 5 is an enlargement of the numerical Schlierens 
around the cavity at time step 225,000. Eddies are 
clearly displayed. 

3.2 Low speed case {v = 26.8m/s) 

In this case, all the conditions are the same as above ex- 
cept the mean flow Mach number is lower, M — .0788, 


corresponding to the car speed of v — 26. 8rn/s. For run- 
ning at such a low Mach number, most of the numerical 
schemes need some form of preconditioning. However, 
this doesn’t seem to bother the CE/SE scheme. Totally 
390, 000 time steps are run for this case and the PSD is 
shown in Fig. 6. The strongest (also the lowest) tone fre- 
quency is estimated to be about 1,911 Hz. Since this 
frequency is not far away from 1 , 824 Hz, while the car 
speed is almost reduced by half, it is guessed that the tone 
may come from the Helmholtz resonant frequency of the 
cavity. 

4 Concluding Remarks 

A typical vortex induced, self-excited tonal gap noise 
is simulated and compared to the experimental data. 
Through the ’difficult’ problem, the capability of the new 
CE/SE scheme is demonstrated. Specifically, 

1. the novel non-reflecting boundary conditions based 
on flux balance is simple, genuinely multi- 
dimensional, and easy to implement, it is effective 
even for near field boundaries; 

2. the (2nd order) CE/SE scheme is robust, efficient 
and yields high resolution, low dispersion results 
similar to those of higher-order schemes. 
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